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ABSTRACT 

We have mapped the superwind/halo region of the nearby starburst galaxy M82 
in the mid-infrared with Spitzer — IRS. The spectral regions covered include the 
fpStl) — S( 3), [Nell], [Neill] emission lines and PAH features. We estimate the total 
warm H 2 mass and the kinetic energy of the outflowing warm molecular gas to be 
between M warm ~ 5 — 17 x 10 6 M 0 and Ek ~ 6 — 20 x 10 53 erg. Using the ratios of 
the 6.2, 7.7 and 11.3 micron PAH features in the IRS spectra, we are able to estimate 
the average size and ionization state of the small grains in the superwind. There are 
large variations in the PAH flux ratios throughout the outflow. The 11.3/7.7 and the 
6.2/7.7 PAH ratios both vary by more than a factor of five across the wind region. 
The Northern part of the wind has a significant population of PAH’s with smaller 
6.2/7.7 ratios than either the starburst disk or the Southern wind, indicating that on 
average, PAH emitters are larger and more ionized. The warm molecular gas to PAH 
flux ratios (H 2 /PAH) are enhanced in the outflow by factors of 10-100 as compared 
to the starburst disk. This enhancement in the H 2 /PAH ratio does not seem to follow 
the ionization of the atomic gas (as measured with the [NeIII]/[NeII] line flux ratio) 
in the outflow. This suggests that much of the warm H 2 in the outflow is excited by 
shocks. The observed H 2 line intensities can be reproduced with low velocity shocks 
(v < 40 km s _1 ) driven into moderately dense molecular gas (10 2 < nn < 10 4 cm -3 ) 
entrained in the outflow. 
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1 INTRODUCTION 

Galactic outflows, “superwinds”, are driven by hot gas 
and momentum generated by the combined mechanical 
and radiative energy and m omentum outpu t from stellar 
winds and supernovae (e.g. I Heckman et al.l 1 1990h . They 
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have been invoked as a source for the heating and metal- 
enrichme nt of both the intra- c luster and inter-g alactic 
medium (lAdelberger et ahl 120031 : iLoewensteinl 12004 ), as a 
major driver of galaxy evolution (e.g. Croton et al. 120061 ). 
and as playing a crucial role in establishing the mass - 
metallicity relation among galaxies (ITremonti et all 120041 1. 
Superwinds have also been hypothesized to be the mech¬ 
anism through which a dust enshrouded ultraluminous in- 
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frared galax y (U LIR G) transitio ns to an optically identified 


irarea galax y ( u mnt j) tra nsitior 
quasar (e.g. ISanders et aPil988l 'l. 


The dynamical evolution of a starburst- 
driven outflow has been extensively d is cussed (e.g . 
Chevalier et al. 1985: Suchkov et al. 1994; Wang 1995; 


Tenorio-Taglc & Munoz-Tunon 1998; Stric kland et al.1 


200df . Superwinds are generated when the kinetic energy 
in the outflow from massive stars and supernovae is 
thermalized, generating a region of very hot (T ~ 10 8 
K] low-density gas in the ISM of a starburst galaxy 
dChevalier et al.lll985l l. This hot gas will expand, sweeping- 
up and entraining ambient gas and creating a hot cavity 
within the overall structure of the ISM. This region of hot 
gas is highly over-pressurized relative to its surrounding 
ambient ISM and will expand along the steepest pressure 
gradients. As it expands, it sweeps up and sh ocks ambient 
mate ri al, creating an outf lowing “superwind” (ICastor et al.l 
1 19751 : IWeaver et al.l I1977T ). Contributing to this overall 
energy driven outflow is the radiation pressure on dust 
grains t hat are lik ely coupled to the ionized gas within the 
flow dMurrav et al.l[2005l'l . 

There is a significant amount of morphological, phys¬ 
ical and kinematic evidence for the existence of super- 
win ds in nearby sta r burs t and infrared luminous galaxies 
(e.g. lHeckman et al.ll 1990l i. Large- scale optical emission-line 
and associated X-ray nebulae are ubiquitous in starbursts 
( [Heckman et al.lll990l : iLehnert et al.lfl999l : IStrickland et al.l 
l2004fh and these increase in size and luminosity from 
dwarf starbu rsts to ULIRGs, w here they can be tens of 
kpc in size dCrimes et all l2()0ol ). In additio n, the winds 
are observed in molecula r lines such as CO (Walt er"et al.l 
120021 : iBolatto et al.l 120131 ') and blue-shifted absorption fea¬ 
tures (e ; g ; _J^aD}_seen_ 2 K)jected_againsWthe starburst nu- 
clei dHeckman et al.l 120001 : iRupke fc Veilleuxl 120051 : iMartinl 
l2005l l indicating that the neutral, cold ISM is being swept- 
up in the outflow. Recently, vigorous outflows have also been 
detected in ULIRGs and LIRGs in molecular absorption 


extra-planar warm H 2 knots and filaments extending more 
than ~ 5 kpc above and below the galactic plane of M82 
have bee n imaged in th e near -infrared using narrow-band 
filters by IVeilleux et al.l (l2009l h Using wide-field imaging, 
they found that H 2 is widespread, suppressed relative to 
PAHs in the fainter regions, and concluded that H 2 is not 
a dy namically important component of the outflow. How¬ 
ever, IVeilleux et al.i J2009l l could only measure the very hot 
1000 K) H 2 gas, and the mass of very hot H 2 gas is much 
smaller than the mass of the H 2 gas probed by this study, 
as we show in Section 3.3. 

In this paper, we study the M82 outflow with Spitzer — 
IRS. With Spitzer/IRS spectroscopy, we can study in de¬ 
tail all the PAH features, derive the grain properties through 
analysis of the PAH ratios, and compare these to the warm 
(200K) molecular gas and the ionized gas (via the strong Ne 
lines) as a function of position above and below the plane of 
M82 for the first time. The rotational lines of H 2 , observable 
in the infrared, unlike the near-infrared ro-vibrational lines, 
probe most of the warm molecular gas mass. We then use the 
ratio of pure rotational H 2 lines to constrain the physical and 
excitation conditions of the gas. We use PAH 11.3/7.7 and 
6.2/7.7 ratio maps to study PAH sizes and ionization across 
the outflow, [NeIII]/[NeII] to study the variation of the hard¬ 
ness of the ionization field, H 2 /PAH to diagnose the origin 
of the H 2 emission, and the H 2 /[NeII] and PAH/[NeII] line 
flux ratios to help understand the survivability of the PAH 
grains in the wind. For the first time we can make robust 
maps of these ratios throughout the outflow of M82, and 
compare all of these to the spatially-resolved X-ray emis¬ 
sion to understand the properties of the dust and molecular 
gas and ionized gas in the wind. 

The paper is organized as follows: in § [2] we describe 
the observations and the data reduction; in § [3] we present 
and analyze the the spectral maps built from the observa¬ 
tions, using the PAH 11.3/7.7 and 6.2/7.7, [NeIII]/[NeII], 
H 2 /PAH, H 2 /[NeII], and PAH/[NeII] line ratios, and also 


lines with Herschel 

— PACS which sug 

gest very high mass H 2 excitation diagrams to diagnose the properties of the ion- 

outflow rates (e.g. 

Sturm et al. 

201h; 

Fischer et al.1 120101; ized gas, PAHs and the H 2 emission; in §|J]we diagnose the 

Veilleux et al. 20131; 

Spoon et al. 

20ll|). ULIRGs which have origin of H 2 emission using X-ray images and shock excita- 


higher AGN bolometric fractions appear to have higher ter¬ 
minal velocities and shorter gas depletion timescales than 
those with lower AGN bolometric fractions, suggesting that 
at least in ULIRGs, the AGN pow er is related to the 
feedback on the c old molecular gas (IVeilleux et al.l 120131 : 
ICicone et alJl2014h . 

M82 is the brightes t, nearest, and best-studied starburst 
galaxy with an outflo w (Heck man et al.lll99d : IShopbell et al.l 
ll99^jLehnert et_al.l Il999l '). At a distance of about 3.5 Mpc 
( Jacobs et al.l 2009l i , M82 affords an unparalleled oppor¬ 
tunity to study a dusty outflow in great detail. The out¬ 
flow in M82 has extensive and complex extra-planar poly¬ 
cyclic aromatic hydrocarbon (PAH) and H 2 emission, ev¬ 
ident from the Spitz er— IRAC and S p itzer — IRS obser¬ 
vations presented in lEngelbracht et aD (l2006l h Extensive 
PAH emission is also seen at 3.3/xm in the AKARI data 
lYamagishi et al] (|20ld h PAHs are thought to be the ori¬ 
gin of strong, broad emission features arisin g in the mid- 
infra r ed spectra of star fo rming galaxies (e.g. iBrandl et al.l 
120061 : ISmith et al.l l2007al l. M82 has also been observed in 
the [CII] and [01] emission lines by IContursi et al.l (120131 1 
who find a low velocity outflow in the atomic gas. Complex 


tion models; and in §[5]we present the conclusions. 


2 OBSERVATIONS AND DATA REDUCTION 

The outflow of M82 was observed with Spitzer-IRS between 
November 10-14 2006 and May 9-22 2008, using all low- 
resolution modules. In Figure[T]we show the footprints of the 
IRS observations on the outflow of M82. We extracted only 
the maps that contain both modules SL1 (7.8 — 14/rm) and 
SL2 (5.5 — 7.8/rm) or LL1 (20 — 35/xm) and LL2 (14 —20/rm). 
In total we mapped an area of approximately 5 aremin 2 

above and below the plan e of M82. _ 

We used CUBISM (ISmith et al.ll2007al ') to build SL and 
LL spectral cubes from from the Spitzer-IRS Basic Calibra¬ 
tion Data (BCDs). In order to create and subtract a local 
background from the wind spectra, we observed a region 20 
kpc North of the M82 disk with both IRS modules. This 
region was observed with the same parameters (integration 
time per pixel and cycles per position) that were used for the 
wind. After verifying that no line, PAH or dust continuum 
was visible in the off field, we subtracted an average image 
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Figure 1. Overlay of the footprints of SL (5 — 14/im - magenta) 
and LL ( 14 — 35/mi - cya n) maps on a Spitzer —IRAC 8^xm image 
dEngelbracht et aD [2QQ6'). The “center” region of M82 mapped in 
iBeirao et alJ (|20Q 3) is presented in green. The scale is 1 kpc per 
arcminute. 

created from this region from the wind data at each location 
and for each IRS module before building th e spectral cubes. 

We used PAHFIT llSmith et al.ll2007bl l to create maps 
of spectral lines and PAH features. PAHFIT is a spectral 
fitting routine that decomposes IRS low-resolution spectra 
into broad PAH features, unresolved line emission, and con¬ 
tinuum emission from dust grains. PAHFIT allows for the 
deblending of overlapping features, in particular the PAH 
emission, silicate absorption and fine-structure atomic, and 
warm H 2 emission lines. Because it performs a simultaneous 
fit to the emission features and the underlying continuum, 
PAHFIT works best if it can fit combined SL and LL spectra 
across the full 5 — 40/rm range. In order to provide a stable 
long wavelength continuum from which to estimate the true 
PAH emission, including the broad feature wings, we formed 
a combined SL + LL cube. In order to investigate the spatial 
variations of the H 2 /PAH ratio at the maximum resolution 
we decided to rotate and regrid the LL cubes to match the 
orientation and pixel size of the SL cubes (l".85 x l".85). 
Thus we provide a stable long-term extrapolation of the LL 
fit to the SL fit. By calculating the flux in the new LL cubes 
using a bilinear interpolation (in surface brightness units), 
we assure that the flux over the native SL pixel scale is con¬ 
served. 

Finally, we applied a small correction factor to each in¬ 
dividual SL spectrum, necessary to match the continuum in 
the corresponding LL spectrum. This scale factor ranges be¬ 
tween 0.65 — 0.99, and was calculated by comparing pixels in 
the spectral overlap region between SL and LL at approx¬ 
imately 14pm. This scale factor is applied before running 
PAHFIT on the combined SL + LL cubes. Such scaling is 
commonly observed in spatially extended sources as a result 
of the different entrance slit dimensions of the SL and LL 



Wavelength 

Figure 2. Detailed PAHFIT spectral decomposition of a 5 — 
20pm spectrum of the M82 outflow. The red solid lines represent 
the thermal dust continuum components, and the thick gray line 
shows the total (dust+stellar) continuum. Blue peaks set above 
the total continuum are PAH features, while the violet peaks are 
unresolved atomic and molecular spectral lines. A fully mixed 
dust extinction model is used (shape indicated by the dotted black 
line at the top), with the relative extinction given on the axis at 
right. The solid green line is the full fitted model, overlaid on the 
observed flux intensities and uncertainties. 

spectrograph modules. In Figure [2] we present an example 
of the spectral decomposition performed at every location of 
the map with PAHFIT. The output of PAHFIT is saved for 
each feature, and is used to create all the maps (PAH, H 2 , 
[Nell] and [Neill]) at the native resolution of the SL data 
as shown in Figure [3] Although care was taken to avoid the 
FIR-bright nucleus of M82 during the spectral mapping, the 
LL1 data exhibit a “ripple” pattern that varies spatially and 
spectrally, and the LL2 data have a slightly elevated contin¬ 
uum level in some pixels (compared to SL1). We believe both 
of these to be the effect of light scattered into the nuclear 
IRS entrance slits. Therefore we have opted not to use the 
LL1 data in the fitting and analysis, except for the deriva¬ 
tion of an upper limit on the H 2 S(0) line flux, calculated 
as the 3 — a noise level at 28/.tm. We have excluded regions 
of the LL2 data cubes where the implied scale factor in the 
continuum between SL1 and LL2 is greater than 2.0. 

3 RESULTS 

3.1 Morphology of the dust, ionized gas and 
warm H 2 emission 

In Figure [3] we present the total PAH emission map calcu¬ 
lated by summing over all the PAH bands between 5 — 14pm, 
the [NeII]12.8/.tm flux map derived from the SL data cubes, 
the total H 2 map derived from the SL and LL2 cubes, and 
the [NeIII]15.6/.tm emission map. We find widespread emis¬ 
sion with significant structure in all of the line maps. To the 
North, both PAH and [Nell] exhibit two prominent, broad 
radial filaments. To the South, the PAH map shows a com¬ 
plex structure with a bubble, or loop. This loop is also seen 
in the [Nell] map. The overall morphologies of the warm H 2 
and Ha are generally very similar, as seen in the NIR by 
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Figure 3. Maps of the continuum-subtracted total PAH emission between 5 — lAfim (upper left), the [NeII]12.8p,m emission (upper 
right), total H 2 S(l)-S(3) emission (lower left) and total [NeIII]15.6^m (lower right). The maps are clipped at 2-sigma level in all cases. 
The white radial strip in the bottom two panels are regions that have been masked out due to scattered light in the LL2 module, which 
affects the H2S(1) and [NeIII]15.6/rm line maps. 


IVeilleux et alJ (|2009|). Maps of t he col d molecular gas made 
in CO (J=l-0) by Salak et all <| 20131 ) are at lower spatia 
resolution, but also show extended, extraplanar gas. 

We can study the physical properties of the ionized gas 
and PAHs by mapping the ratios between ionized lines such 
as [NeIII]/[NeII] and between the different PAH bands. If 
the ionized gas is photoionized, the Ne line ratio is sensi¬ 
tive to the hardness of the radiation field and the ionization 
parameter. The Ne ratios observed in the outflow are stan¬ 
dard values for the hardness of the radiation field found in 


starburst galaxies. The ratio of the 6.2, 7.7 and 11.3 micron 
PAH features can be used to estimate the ionization stat e 
and the size of the small dust grains fe.g. iDraine fe Tl200ll') . 
Even the overall absence of PAH emission may be related 
to the ionization state of the gas since we do not expect gas 
emitting strongly in [Neill] to show much PAH emission, as 
the small grains are easily destroyed by energet i c photons 
(e.p. lEngelbracht et al . 1 1200 ( 3 ; iMadden et al .1 1200 ( 3 :IWu et all 

l2006t l. The 11,3//m/7.7//m ratio is high when PAHs are more 
neutral and larger, while the 6.2pm/7.7/rm PAH ratio is high 
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Figure 4. Maps of the PAH 11.3/7.7 ratio (upper left), PAH 6.2/7.7 (upper right), H2[S(1) — S(3)]/PAH(7.7 + 8.6pm) (lower left) and 
[Neill]/[Nell] (lower right). All maps have been clipped at the 3 — a surface brightness level. The white radial strip in the bottom two 
panels are regions that have been masked out due to scattered light in the LL2 module, which affects the If 2 ,S’( 1) and [Neill] lines. The 
circles represent the selected regions where the spectra in § [5] were extracted and analyzed. The fluxes are presented in Table 1 and 2 
(see §2 for details). 


when PAH grains are smaller. However, we caution that vari¬ 
ations in the PAH line ratios are also influenced by the struc¬ 
tural properties and elemental abun dances in small particle s 
giving rise to the PAH features (e.g. lYamagishi et al.ll2012l) . 

In Figure [4] we present maps of the PAH 11.3/7.7 and 
6.2/7.7 ratios, the H 2 /PAH ratio and the [NeIII]/[NeII] ra¬ 
tio. The H 2 /PAH ratio map was produced using a H 2 [S(l)- 
S(3)] total flux map and PAH (7.7 + 8.6pm) flux map. Al¬ 
though there is a great deal of structure in the PAH ratio 


maps, the general trend is that the 11.3/7.7 and the 6.2/7.7 
ratios increase as we move away from the starburst disk. The 
PAH 11.3/7.7 is enhanced by a factor of ~ 4 and the PAH 
6.2/7.7 is enhanced by a factor of ~ 3 at the northern edge 
of the outflow compared to the regions near the center of the 
starburst disk. There is an enhancement of the 11.3/7.7 and 
6.2/7.7 ratios by ~ 2 in the southeast edges of the outflow. 

The H 2 /PAH ratio varies between 0.01 — 1, and the 
[NeIII]/[NeII] ratio varies between 0.2 — 0.9. In the region 
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between the PAH emission “fingers” (as particularly evident 
in the Northen outflow region in Figure^, both maps show 
an enhancement of the ratio. 

In Figure [4] we also mark the 12 regions from which 
we have extracted 5 — 20pm spectra and performed multi¬ 
temperature fits to the H 2 lines. We show these spectra in 
Figure [5] These regions were chosen to represent the full 
range of H 2 /PAH ratios and to approximately sample the 
full spatial extent of the outflow. For comparison, we also 
sh ow an integrat ed spectrum of the center of M82 (Figure(5] 
IBeirao et aljfeoosh . which we use as a reference. In addition 
to the 12 regions, we also show a composite spectrum of 
the entire IRS map in Figure [5j We measure the fluxes of 
the PAH features between 5 — 14pm, the fine-structure lines 
[Nell] and [Neill], and the rotational excitation lines of H 2 
from S(0) to S(3) and present the results in Tables 1 and 2. 
Regions 1, 5, 6, 7 and 9 are the most similar to the reference 
spectrum in slope and the relative strength of PAH features, 
which is not surprising since they are located close to the 
center. PAHs are also strong in regions 3, 8, 10, 11 and 
12, but their spectrum also exhibits a flatter continuum, 
probably indicating the presence of colder dust. Region 3 is 
the most distant from the central starburst and it has the 
weakest PAH emission. Region 4 exhibits a flat continuum, 
but a comparatively strong PAH 11.3pm feature. 

As we have seen in the maps, the spectra show strong 
variations of the H 2 /PAH ratios between the regions. The 
H2S(1) line is very prominent in the spectrum of Region 3 
compared to the PAH bands, while in the regions close to 
the disk, such as Regions 7 and 9, this line is much fainter 
compared to the PAH features. Indeed, using the values 
in Tables 1 and 2, we calculate in Region 3 a very high 
H 2 /PAH ~ 0.08, while in Region 9, H 2 /PAH ~ 0.01. 

I 11 Figure[6]we present a plot of the PAH 11.3/7.7 ratio 
as a function of the 6.2/7.7 PAH ratio. The range of 6.2/7.7 
PAH ratios in the outflow is much larger and the values 
are, on average, smaller than the disk. The outflow exhibits 
6.2/7.7 and 11.3/7.7 ratios ranging between ~ 0.2 — 0.4. 
On average, the bulk of the superwind shows larger, more 
ionized grains than in the starburst disk, although there are 
regions where the ratios overlap. As an example of the differ¬ 
ent PAH ratios found in the outflow, we present in the upper 
right corner of Figure [6] two spectra from different regions of 
the outflow, and an average spectrum of the disk. In general, 
the outflow has stronger 17pm PAH and H 2 emission than 
the SB disk, and the spectrum in the region 2 kpc north of 
the disk (the blue spectrum shown in Figure O has higher 
17pm feature and H 2 S(l) flux, indicating the prevalence of 
larger PAH grains and enhanced H 2 emission further away 
from the disk. In summary, we see a much larger range in 
the PAH ratios in the outflow than is seen in the starburst 
disk. The average spectrum of the wind has smaller 6.2/7.7 
ratios and larger 11.3/7.7 ratios than the average spectrum 
of the starburst, suggesting the grains may be larger and 
more ionized, on average, than those in the starburst disk. 


3.2 H 2 /PAH variations in the M82 wind 

In Figure|7]we plot the variation of H2S'(1) — S(3)/PAH7.7+ 
8.6pm as a function of [Nell] surface brightness. Since the 
[Nell] emission arises predominantly from the starburst disk, 
this traces the relative H 2 and PAH emission as a function 




Figure 5. Mid-infrared 5 — 20pm spectra from selected regions 
i n the outflow, as seen in Figure [3] and the center of M82 from 
IBeirao et alj ll2008f) . The spectrum labeled “Total” is the sum of 
all pixels in the map. 


of distance fr om the starburst. This ratio has been previ¬ 
ously used by I Ogle et all (1201011 and iGuillard et al.l (I2012fl 
to show the relative importance of the UV versus mechani¬ 
cal heating of the gas. We can see that the H 2 /PAH ratio 
decreases with increasing [Nell] flux. The vast majority of 
the mapped regions lie well above the values seen in star- 
forming galax i es in t he local Universe (the dashed lines - see 
iRoussel et al.l (120071 41. The H 2 /PAH ratio values span two 
orders-of-magnitude, 0.01-1. Most of the M82 wind points 
have H 2 /PAH between 0.01 and 0.10, but ~ 4% have ex¬ 
treme values, H 2 /PAH > 0.5. The fact that the H 2 /PAH 
ratios seen in the wind are well abov e those seen in the SB 
disk, the normal gala xies observed b y Roussel et al. ] teOOTf) . 
and PDR models of IGuillard et al.l (120121 4 (see section 4) 
suggests that the H 2 is not solely excited by stellar radia¬ 
tion from the starburst. 

We can also investigate possible correlations of the 
H 2 /PAH ratio with diagnostics of the atomic gas and the 
size and ionization state of the grains. In Figure [8] we show 
H 2 /PAH as a function of the PAH 11.3/7.7 ratio, the PAH 
6.2/7.7 ratio, the H 2 S(3)/S(l) ratio, and the [Neill]/[Nell] 
ratio. We see that in the disk the H 2 /PAH ratio increases 
with increasing PAH 11.3/7.7 ratio. While the trend of in¬ 
creasing ^ 2 /PAH in the wind is similar to that seen in the 
disk, the wind values extend to much larger values of both 
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Figure 6. Plot of the PAH ratio 11.3/7.7 as a function of 6.2/7.7 
for 10 // x 10 ,/ regions in the outflow. Red represents the north¬ 
ern outflow, cyan the sout hern outflow and g reen represents the 
central 1 kpc in M82 from IBeirao et alj <120081). The dashed lines 
represent neutral and ionized PAHs from lDraine &; O (|200ll) . The 
dot-dashed line represent the PAH ratio values for a PAH grain of 
250 carbon atoms. The red and cyan crosses represent typical un¬ 
certainties in the north and southern outflows. The three spectra 
on the upper right corner are derived from different regions in the 
outflow (in units of MJy/sr). The blue spectrum is derived from 
a region 2 kpc north of the disk and is represented by the blue 
data points, the red spectrum is an average spectrum of the whole 
northern outflow, and the green spectrum the average spectrum 
of the central region. 



11.3/7.7 PAH and H 2 /PAH than seen in the starburst. The 
scatter to high H 2 /PAH and a much larger range in the 
6.2/7.7 ratio is seen in the upper right panel of Figure [8] 
as well. In the lower left, we observe an anti-correlation be¬ 
tween H 2 /PAH and S(3)/S(l) in the disk as warmer gas 
has lower H 2 /PAH. The wind points seem consistent with 
this trend but the uncertainties are large. This general lack 
of very warm molecular gas and high H 2 / PAH ratios in the 
wind may be a function of inability of the MIR diagnostic 
features to probe the hottest gas. The [Nelli]/[Nell] also in¬ 
creases with H 2 /PAH in the disk, and the wind points to 
the North seem consistent with the trend, but clearly there 
are many wind points in a few regions in the northern out¬ 
flow which have much higher H 2 /PAH ratios than seen in 
the disk at a given [NeIII]/[NeII] value. 

To further help us understand what drives the vari¬ 
ations of the 112/PAH flux ratios, we present maps of 
H 2 /[NeII] and PAH/[Nell] in Figure [9] In general, the re¬ 
gions with high H 2 /PAH correspond to high PAH/[Nell], 
suggesting that it is not simply a decreasing PAH that causes 
the rise in H 2 /PAH off the plane of M82. Only in the north¬ 
ern edge of the map do we see a lower PAH/[Nell] ratio, 
which might indicate reduced PAH heating efficiency. It is 
important to note, however, that while the PAH/[Nell] ratio 
does change by factors of a few across the map, this is much 
less than the variation observed in H 2 /PAH. Therefore the 
latter cannot be solely explained by a change in the strength 
of the PAH emission alone. 


1.000 
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Figure 7. Plot of the H 2 S(1)-S(3)/PAH (7.7 + 8.6/im) as a func¬ 
tion of the [Nell] surface brightness. Each datapoint represents 
a 10 // x 10 // region. Red points are the northern outflow, cyan 
the southern outflow, and green the central 1 kpc of M82 from 
IBeirao et alj <l2008l) . The dashed line and the dash-dotted lines 
represent the average and the la variation, respectively, of the 
H 2 S(1)- S(3)/PAH (7.7 + 8.6/i m) ratio in the star-forming galaxy 
sample of lRoussel et al.l <|2007|) . The dotted-dashed line represents 
the upper limit in the H 2 S(1)-S(3)/PAH (7. 7 + 8.6pm) ratio pre¬ 
dicted by PDR models (iGuillard et alj|2012h . 

The error bar in red represent typical uncertainties in the north¬ 
ern outflow, while the typical uncertainties in the data for the 
southern outflow region are smaller than the points. 


3.3 H 2 masses and temperatures 

With the rotational H 2 emission line fluxes in Table [2[ we 
can estimate the gas temperature, column density, and mass 
distributions, assuming that the ortho/para ratio is in local 
thermal equilibrium (LTE). In Figure ITOl we plot the exci¬ 
tation diagrams for all regions, except 2 and 9, which did 
not have sufficient signal-to-noise for fitting the data with 
these sorts of models. Using the Boltzmann equation, we 
fit the H 2 lines with two components, each with different 
temperatures. We then derive the total column density, H 2 
temperatures and mass distributions. In fitting the data we 
set the minimum allowed temperature to 100 K and treated 
every upper limit as a 3-sigma detection. The S(l) line is the 
line with the highest signal-to-noise. Therefore, when we re¬ 
move the S(0) line as a detection, we obtain essentially the 
same results. 

The diagrams for Regions 1, 3, 4, 5, 7 and 11 look very 
similar in shape, derived from very similar temperatures for 
each component. In Region 8, both components have similar 
temperatures, while in Region 10, there is a large difference 
in temperature between the two components. We note that 
Regions 8 and 10 have the most discrepant S(0) upper limits 
relative to the other fits, making the mass estimates partic¬ 
ularly uncertain. In Tabic [3] we present the properties of the 
warm H 2 gas in these regions of interest. In all the regions 
the temperature of the warmer component is above 500 K, 
especially w arm when compare d with typical star forming 
regions fe.g. lRoussel et aOl2007l 'l. The column density of the 
bulk of the warm H 2 gas ranges by a factor of ~ 30. The 
mass ratio between the warmer and colder component varies 
between 0.0003 — 0.02, with the lowest value observed in Re¬ 
gion 8 and the highest in Region 1. The region with the most 
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PAH 11.3/7.7 


PAH 6.2/7.7 




Figure 8. Plot of the H 2 S(3)-S(1)/PAH (7.7 + 8.6^m) as a function of PAH 11.3/7.7 (upper left), PAH 6.2/7.7 (upper right), H 2 
S(3)/S(l) (lower left) and [Neill]/[Nell] (lower right). Each datapoint represents a regi on of 10 /; x 10 /; in extent. Red points are the 
northern outflow, cyan the southern outflow, and green the central 1 kpc of M82 from [Beirao et all i 20081) . The error bars represent 
typical uncertainties in the northern (red) and southern (cyan) outflows. The d ashed lines represent the average and the lcr variation of 
the H 2 S(1)-S(3)/PAH (7.7 + 8.6/im) ratio in the star-forming galaxy sample of [Roussel et al. (2007)- T he dot ted-dashed line represents 
the upper limit in the H 2 S(1)-S(3)/PAH (7.7 + 8.6/rm) ratio predicted by PDR models dGuillard et al.l[2012t) . 


luminous H 2 emission is Region 8, where the warmer compo¬ 
nent accounts for ~ 95% of the luminosity, whereas the next 
most luminous region, 11, the colder component accounts for 
~ 65% of the luminosity. 

We also integrated the emission over the entire observed 
nebula surrounding M82 and measured the line fluxes, which 
are included in Table [2] The total warm H 2 mass derived 
from the integrated emissio n is ~ 4.8 x 10 6 Mg. For com¬ 
pleteness, we note that IVeilleux et al.l d2009i ') finds a mass 
of H 2 at a temperature T>1000 K (thus hotter than our 
“hot” gas) of 1.2xl0 4 Mg, based on near-infrared obser¬ 
vations. The relatively low mass is not surprising since the 
near-infrared ro-vibrational lines are probing very high ex¬ 
citation gas which emits very efficiently in these lines. The 
mass of warm molecular gas in the wind that we are detect¬ 
ing with the IRS is thus more than a factor of 100 more 
than is estimated from the NIR line maps. We have also es¬ 
timated the H 2 masses from fitting the H 2 line fluxes with 
molecular shock models (see Section 4.3). 

Considering a cold molecular ga s mass of ~ 10 9 Mq 
( Walter et al. ' ’20021; ISalak et ali 12013 ') , we calculate a ra¬ 
tio of warm to cold molecular gas in the outflow of about 
0.01. Therefore the warm molecular gas we measure here is 


a relatively minor component of the total molecular mass in 
the ISM. Care must be taken when making this comparison 
as the emission regions for the CO and H 2 do not overlap 
completely. In fact, the most extended CO emission is only 
within a few arc minutes of the plane a nd much less extended 
llSalak et al.ll20l3 : IWalter et al Jl20()2h , although new studies 
trace the CO-emitting gas to larger distances dLerov et al.1 
20 i li) . However, we note that most of the estimated mass 
for the warm molecular gas in our observations also comes 
from regions which lie close to the disk (and likely the base 
of the wind). 


4 DISCUSSION 

4.1 H 2 heating processes 

In the scenario of thermal excitation of H 2 , three mecha¬ 
nisms dominate the heating: (1) UV radiation from star- 
burst, (2) X-rays from the starburst and win d, and /or (3 ) 
shocks induced by the outflow. Following lOele et al . (2010 '). 
we use the H 2 -to-PAH luminosity ratio to investigate the 
contribution of UV photons to the total heating of the H 2 














































Spitzer-IRS Maps of the Superwind in M82 9 



Right Ascension Right Ascension 


Figure 9. Maps of the H2>S(1) — S'(3)/[Nell] ratio (left), PAH5 — 17/rm/[NeII] (right). H 2 is the sum of the S(l), S(2), and S(3) H 2 line 
fluxes and PAH5 — 17//rn is the sum of the fluxes of all PAH features between ~ 5 and 17//rn. All maps have been clipped at the 3 — rr 
surface brightness level. The white radial strip in the bottom two panels are regions that have been masked out due to scattered light in 
the LL2 module, as in Figures [3] and [4] 


gas in the outflow. In starburst galaxies like M82, we expect 
the H 2 emission to be predominantly powered by UV radi¬ 
ation from young stars but undoubtedly, processes such as 
low-velocity shocks and turbulent dissipation contribute to 
the H 2 emission. It should also be noted that II 2 /PAH ra¬ 
tios on the order of 0.01 have also been seen from translucent 
clouds in the cold, diffuse ISM of the Milky Wa y, w hich are 
thought to be excited by UV photons (jlngalls et al.ll201l[ ). 

There is a clear offset of the wind points from those 
in the SB disk in terms of H 2 , PAH and ionized atomic 
emission (Figs. 7, 8, & 9). The data in regions of the out¬ 
flowing gas show a large range in both [NeIII]/[NeII] and 
H 2 /PAH, but are clearly offset to higher H 2 /PAH at a 
given [NeIII]/[NeII] ratio compared to regions within the 
SB disk, suggesting an additional excitation mechanism at 
work. In fact, models suggest that the large values of the 
H 2 /PAH ratios observed in so me regions (>0 .1) cannot be 
explained by photo-ionization (IGuillard et al.ll2012l ). 

The H 2 gas could be indirectly excited by X-rays 
emitted from the expanding bubble of hot and tenuous 
;as generat ed by the intense star formation within M82 
Begelman et al. 119891 ). high velocity shocks in the en¬ 
trained gas, X-rays due to inverse Compton scattering of 
IR photons. In M82, the X-ray emission covers approxi¬ 
mately the same region a s the optical emissio n line fila¬ 
ments llWatson et al.fl 19841: IShonbell et, al.lll998l ) and radio 
emission ( Seaauist fe Odegard Il99ll ). For example, the high 
H 2 /PAH in the northern tip ~5 kpc from the disk, appears 
to be associated with a region with X-ray emission (Fig¬ 
ure [TT]). 

We can use the total H 2 5(l) — 5(3) luminosity and the 
total X-ray luminosity in the outfl ow to c alcu late the H 2 - 
to-X-ray luminosity ratio (IStevens et al] 120031 ). The H to 


H 2 transition can occur at relatively low column densities, 
N h «10 21 cm -2 . At such low columns, even relatively soft 
X-rays can penetrate. Our estimates for the M 82 outflow are 
L(H 2 S( 0) - 5(3)) /L x (2-10 keV) - 1 (IStevens et al.ll2003D 
and L(H 2 S(0) - 5(3)) /L x (0 2 - 4k eV) ~ 0.5 llWatson et al.l 
11984 IStrickland fe Heckmanll2007l) . over the northern and 
south ern IRS map. Using XDR models (e.g. iMalonev et al.l 
Il996h we can estimate whether the H 2 line emission could 
be powere d by X -ray heating, following the procedure de¬ 
scribed in I Ogle et alJ (l20ldl . In these models, 30% — 40% 
of the absorbed X-ray flux goes into gas heating via pho¬ 
toelectrons, and the cooling by H 2 rotational lines in XDR 
models is ~ 2% of the total gas cooling for a gas temperature 
of 200 K. At this temperature, the ratio of the first four ro¬ 
tational lines to the total rotational line luminosity is L(H 2 
0-0 S(0)-S(3))/ L(H 2 ) = 0.6. Combining the above factors, 
we estimate a maximum H 2 to X-ray lumin osity ratio o f L 
(H 2 O-O S(0)-S(3))/ L x (2-10 keV) < 0.01 llGuillard et al.l 
l2012bl ). This upper limit on the H 2 -to-X-ray luminosity is 
conservative because it assumes that all the X-ray flux is 
absorbed by the molecular gas, which is likely not the case 
(since the medium in which X-rays propagate is inhomoge¬ 
neous). Therefore, we conclude that X-rays are not the dom¬ 
inant source of heating of the H 2 gas in the M82 wind, since 
the observed ratios are well above those found in XDR mod¬ 
els. This is also seen in some extreme environments, such as 
powerful radio galaxies, where H 2 /PAH ~ 0.5 — 1 and the 
H 2 /X- ray ratio al so exceeds that expected from XDRs alone 
dOgle et al.ll2010l ). 

Ambient gas clouds heated by the mechanical energy of 
the wind and dense shell fragments from the wind bubble 
carried along by the w ind are the two m ain sources of X-ray 
emission in the wind (IStrickland et al1 E997). In Figure [lT] 
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Figure 10. H 2 excitation diagrams of most selected regions. The dotted lines show a 2-temperature LTE fit of the data. The upper 
limits of the H 2 lines represent a 3-cr noise level at the corresponding line. 




















Spitzer-IRS Maps of the Superwind in M82 11 



09 h 55 m 20 s 10 s 55 m 00%4 m 50 s 40 s 30 s 20 s 


Right Ascension 



09 h 55 m 20 s 10 s 55 m 00%4 m 50 s 40 s 30 s 20 s 

Riqht Ascension loq i/F 

-1 

W m ! sr 


Figure 11 . Image of the H 2 /PAH (top) and total PAH flux 
(5.2 —17pm) (bottom), overlaid with contours from soft (0.3 —1.6 
keV) X-ray emission map of M82 (Strickland & Heckman 2009). 

we show the contours of soft X-ray emission from the outflow 
of M82 overlaid on the H 2 Z.P 2 l.ff and total (5 — 17pm) PAH 
maps. The X-ray and PAH maps show a remarkable cor¬ 
relation, particularly in the southern part of the bi-conical 
outflow, although both maps are projections of the true 3D 
structure and the precise locations of the hot gas and dust 
grains are unknown. Moreover, there are regions of high sur¬ 
face brightness X-ray emission with relatively little PAH 
emission (in the north along the outer wind axis). We also 
see that the middle lobe of X-ray emission in the northern 
outflow, at about 5 kpc from the disk, appears related to 


the increase of H 2 /PAH 7.7 + 8.6pm surface brightness. The 
morphology of H 2 /PAff closely resembles the morph ology 
of the H 2 (2.12pm)/PAH maps in lVeilleux et, al.1 ([2009) . The 
regions where H 2 /PAH >0.1 correspond to regions where 
H 2 (2.12pm)/PAH> 0.003. 

While there is a reasonable spatial correlation between 
the surface brightness distribution of the soft X-ray and H 2 
emission, a correlation is apparent between soft X-rays sur¬ 
face brightness and H 2 /PAff ratio distributions only in high 
surface brightness X-ray emitting regions. There are several 
possibilities why this might be so, but the most obvious is 
that both the X-ray and H 2 emission are fundamentally pow¬ 
ered by the same mechanism, namely, the very hot plasma 
that is acting as the “piston” driving the wind and heat¬ 
ing clouds entrained in the wind. However, the warm H 2 
and hot X-ray emitting plasma should not be co-spatial on 
small scales, and instead the H 2 is likely seen in projection 
against the hot wind, along the edges of the cone in dense 
knots of gas, heated by shocks (see below). Such a situa¬ 
tion would also arise in gas with a wide range of densities or 
through changes as individual clouds are shr edded as th ey 
are a ccelerated by the action of the piston (ICooper et al.l 
120081) . This is also shown by the strong correlation between 
the X-rays and the PAH emission (Figure fill )), especially in 
the Southern wind. As the PAHs and H 2 molecules cannot 
be expected to survive within the 10® —10' A' X-ray emitting 
plasma, it is most likely that they lie along the edges of this 
phase of the outflow. So while the power source is likely the 
same, the two emission regions, not surprisingly, represent 
different physical regions and gas states b e ing im pa cted b y 
the driving very hot plasma (IStrickland fc Heckmanll20071) . 

4.2 The velocities and kinetic energy of the 
molecular gas 

Kinematics of the gas can be used to determine if the differ¬ 
ent gas phases observed in the outflowing wind are related. 
If the observed kinematics of the warm molecular gas as 
probed by the H 2 rotational lines is similar to that of the 
cold molecular gas as probed by the CO lines, the IR atomic 
neutral and ionized emission lines, or the lines of the warm 
ionized medium, then this might suggest a direct physical 
rela tionship between t hese phases. 

ISalak et al.l (i2013l ) derived a velocity of 160 km/s for the 
CO(l-O) emitting gas for the northern outflow region. How¬ 
ever, the farthest extent of the CO emission is only about 
1 arc minute north of the disk of M82, while the IRS map 
extend s to 4 arcmin. The velocity derived by ISalak et alJ 
(l2013ll is not very different from the velocities in the neutral 
and ionized atomic gas in the inf rared within the inne r wind 
region < 1 arcmin from the disk (iContursi et aljfeoish . If we 
restrict ourselves to the inner wind region within a few arc 
minutes of the plane of M82, the velocities in the warm ion¬ 
ized gas are only a few 100 km s -1 or less within a c ouple of 
arc minutes of the disk of M82 llShopbell et al.lll998l ). These 
velocities are consistent wi th those observed in IR neutral 
and ionized atomic lines llContursi et al.l 120131 ). However, 
farther out of the plane of M82, much higher outflow veloc¬ 
ities (~6 00-800 km s -1 ) are o bserved in Ha for the outflow 
in M82 (IShopbell et al.l fl998). Line splitting i s observed in 
both the CO and Ha with similar velocities dWalter et al.l 
120021) . making the association between the cold molecular 
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Figure 12. Shock excitation diagrams of most selected regions in the outflow. The best fit models for either J- (light blue) or C-shocks 
(dark blue) or a combination of each are indicated in the legend to each panel (see text for details). The sum of or the individual the 
best fit models is shown as a solid red line. The level populations of the various rotational levels of the (0-0) H 2 lines are indicated by 
plus signs if they are detection or downward pointing arrows for upper limits and are labeled as S(l), S(2), S(3), S(4), and S(5). 
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gas and the warm io nized gas stronger (cf. ISalak et al .11 201. 'j ; 
IShopbell et a,l .1 f~i~998l) . 

Such a direct relationship would be consistent with 
a picture where the emission we observe in the extended 
outflow of M82 is mainly from entrained and swept-up 
clouds that are predominately shock-heated by a high ram- 
pressure, hot plasma generated within the starburst region. 
If we envision the regions of warm H 2 emission as being the 
areas in the clouds that are shock-heated by the expanding 
plasma, we would expect these p hases to co-exist with the 
warm ionized gas llHeckman et al.lll99Cll l and the cold/warm 
molecular gas. 

Assuming that the velocity of the warm H 2 is similar 
to the d eprojected v elo city of the cold molecular gas de¬ 
rived by ISalak et all (|2013l l, we can calculate its total ki¬ 
netic energy as Ek = 0.5 x MH 2 Voutfiow ~ 5.9 x 10 53 erg 
for Mh 2 ~ 4.8 x 10 6 Mq and v 0 utfiow =160 kmjU 1 . Our es¬ 
timates are considerably higher than those of lVeilleux et al.l 
d2009lf because both their mass estimates are a factor of 500 
less and they assumed a velocity of 100 km s -1 for the warm 
molecular emission. However, they are smaller compared to 
those esti mated for t he co ld molecular and warm ionized 
medium. ISalak et al.l (|2013[ ) estimate that the kinetic energy 
of the wi nd i n the cold molecular gas is « 1 — 4xl0 56 erg. 
IWalter et al.1 (120021 ) estimate a kinetic energy of «3xl0 55 
erg, but they do not map the CO distri bution out to l arge 
distances from the disk compared to lSalak et al.l (|2013h . 

4.3 H 2 supersonic turbulent heating and shock 
models 

The hot plasma driving the outflow will likely transfer 
part of its ram and thermal pressure to swept-up and en¬ 
trained dense molecular gas by driving shoc ks a n d by tur¬ 
bulen t mixing between hot and cold gas (see lGuillard et aid 
120091 1. To model such shocks, we use a grid of stationary 
shock models comp uted with an updated version of the 
iLesaffre et aid d2013l l code to derive the physical shock pa¬ 
rameters that explain the H 2 line ratios. Our initial con¬ 
ditions are that the gas is diffuse, magnetized with a pre- 
shock magnetic field of intensity B\fiG\ = \Jnn [cm -3 ], 
and irradiated with a radiation held of Go = 1 in Habing 
units. As the “Habing held” would photodissociate the H 2 , 
we assume some self-shielding. In the models, we integrate 
the extinction along the model from the pre-shock, and we 
take into account the shielding and self-shielding, which is 
paramet rized i n the photo-dissociation reaction rates of H 2 
and CO llLesaffre et al]|2013t l. At a density of 10 4 cm 3 , the 
ionization fraction in the pre-shock gas is 6 x 10 -8 . This 
matches approximately our estimate of the ionization frac¬ 
tion of 10 -7 at 10 kpc from the X-ray source. We ht the H 2 
line fluxes observed in the different regions with one shock 
or a combination of two shock velocities depending on how 
many H 2 lines are detected, for three diffe rent pre-shoc k 
densit ies: uh = 10 2 , 10 3 and 10 4 cm -3 , f see lGuillard et ahl 
(l2009l'l for details of the method). The results of these fits 
are shown in Figure 1121 where we present the H 2 excita¬ 
tion diagrams for the best velocity combinations (lowest 
X 2 ), and for the intermediate pre-shock density of iih = 10 3 
cm -3 . We present in Table[T]the results of the shock models 
for each region and the integrated spectrum (labeled “To¬ 
tal”). The total gas mass derived from the shock models is 


Mh 2 = 1.7 X 10' Mq, which results in an estimate of the ki¬ 
netic energy of Ek ~ 2 x 10 54 erg. This is a factor of 3 higher 
than the estimates calculated in the temperature diagrams 
of Figure [lO] This difference arises due to the fact that we 
are integrating the postshock temperature profile down to 
50 K (where the H 2 mid-IR lines start to emit), and because 
local thermodynamic equilibrium (LTE) is not assumed in 
the shock models. 

For all regions, the best-fit shock velocities are in the 
range 4-35 km s -1 (grid spans 3-40 km s -1 ; Tabl[4])- These 
velocities are much lower th an the estimated bulk outflow 
velocity llLehnert et al .1 [T 99 il l . The regions show a range of 
best fitting models, but typically the best fitting pre-shock 
densities are nH=100-1000 cm -3 . Our modeling is not in¬ 
tended to be exhaustive, but to simply provide an estimate 
of the physical parameters in the warm molecular gas of the 
outflow. 

At the range of preshock densities and magnetic field 
strengths considered, the C-shocks turn into J-shocks when 
the velocities are above ~ 22 km s -1 . For some regions the 
H 2 line fluxes are reproduced by a combination of two C- 
shocks, in other regions the C-shocks (J-shocks) do not suf¬ 
ficiently populate the highest (lowest) energy levels. In those 
cases, we fit the H 2 line fluxes by a combination of C- and J- 
shocks. The shocks travel through the clouds in ~ 10 4 — 10 5 
years (for 1 pc clouds), which is mu ch shorter than the d y- 
namical time of the outflow (10 s yrs - iLehnert et al.l (jl999t H. 
Using the estimates of the warm molecular gas mass and the 
typical cooling time of the shocked H 2 gas ( t coo i ~ 10 4 yr for 
uh = 10 3 cm -3 ), we can estimate the mass flow of shocked 
gas, M war m/t CO oi ~ 7.5 x 10 6 M q / 10 4 yr = 750M© yr -1 , 
assuming that all the warm H 2 is shock-heated. The dy¬ 
namical timescale of the H 2 -rich outflow is ~ 4.5 kpc/200 
km s -1 = 2.5 x 10 7 yr implying a total mass of shocked gas 
~ 5 x 10 9 Mq. This estimate is likely to be an upper limit 
since the H 2 gas may not extend over the whole outflow 
and because this assumes that the gas is only shocked once, 
which is highly unlikely. The high mass flow rate of shocked 
gas may require the reforma t ion o f H 2 in the post-shocked 
gas fsee iGuillard et ahll2009l . l2010l l to prolong the surviva l 
time of the H 2 gas in the outflow (e.g. [Fragile et alj|2004h . 
Clouds may also survive longer if the radiative cooling time 
is shorter than the time scale necessary for the shocks to 
cross the cloud and the g row th rate of Kelvin-Helmholtz in¬ 
stabilities (ICooper et al.l f2008l. Such processes would lower 
the total molecular mass required to explain the molecular 
emission from the outflow, although for the C shocks con¬ 
sidered here and even for many of the J shocks as well, there 
is little destruction of H 2 . At any rate, the total gas mass 
that is shock heated is similar to that of the cold molecular 
gas. Therefore the reservoir for the shocked gas is likely the 
cold molecular gas as probed through the CO emission from 
the wind. 

4.4 PAH destruction 

Since we see that the PAHs grains are on average larger and 
more ionized in the wind, it is natural to ask whether or 
not some of the enhanced H 2 /PAH could be do to a re¬ 
duction in the PAH emission, rather than an enhancement 
of H 2 . The decrease in PAH flux is probably just due to a 
decrease of the radiation field, but there is also a change 
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in the grain properties associated with the wind. The gra¬ 
dient of decrease in H 2 emission is significantly lower than 
that of the PAH emission, suggesting an additional heat¬ 
ing source for the H 2 (see Section 3.2). One of the ways 
PAH emission can decrease relative to H 2 is PAH destruc¬ 
tion. PAH grains can be destroyed either by hard radiation 
field, X-rays, or shocks. The EWs of the PAH features in 
starburst galaxies are observed to decrease sharply when 
[NeIII]15.6/nn/[NeII]12.8^im> 2, where [NeIII]/[NeII] traces 
the hardness of the radiation field. This has been inter¬ 
preted as evidence of PAH destruction bv intense radiation 
fields llWu et al.ll2006l : [Beirao et al.l 120061 1. In M82, despite 
an increase of [Neill] 15.6/rm/[Nell] 12. 8 /rm in the outflow, 
in most of the outflow [NeIII]15.6/im/[NeII]12.8/.tm remains 
less than 1. [NeIII]15.6/im/[NeII]12.8/rm~ 1 indicates that 
there is hard (> 40.96 eV) radiation present, but not enough 
to have an effect on the average PAH properties. Further¬ 
more, only a small fraction of the gas in the M82 outflow 
is photoionized, and PAHs may survive in HI and H 2 gas. 
PAH destruction is also very sensitive to the clumping fac¬ 
tor of the gas. If the gas is dense and clumpy with relatively 
small volume filling factor and large covering fraction, the 
PAH grains will be shielded from the hard photons in the 
radiation field. So PAH destruction by hard photons is not 
likely to be the cause of the increase of H 2 /PAH in the wind. 

PAH grains are rapidly destroyed in hot, X-ray emitting 
gas dMicelotta et alj|2010br ). However, the presence of PAHs 
at large distances (10 kpc from the disk) in the outflow shows 
that they are protected in cooler phases of the gas, mostly in 
molecular gas. It is possible, however, that X-rays from the 
hot outflow fluid could destroy the PAH carriers, but the 
correlation between the X-ray emitting gas and the PAH 
maps, particularly in the Southern outflow, suggests that 
grain destruction by X-rays in the wind is not dominant. 

Thermal collisions with very small and big grains can 
have a significa nt impact on the co oling of the hot gas ejected 
into the halo (Dwck ct al.1 119901 1. However, the impact of 
this cooling will generally be short lived since the cooling 
efficiency drops as th e grain sputtering occurs in the hot gas 
llGuillard et al.ll20091 1. unless the hot halo gas is constantly 
replenished with dusty gas. A physical dust cooling model 
is needed to fully assess whether thermal emission from col- 
lisionally heated dust in the hot plasma is a significant con¬ 
tribution to the gas cooling in the halo of M82, and whether 
this contribution changes as a function of the radius. 

Another possible mechanism for PAH destruction is 
shocks. As seen in Figure 1121 the H 2 emission that we see 
can be fit with shock models having speeds of up to 40 km 
s . PAH destruction like ly only becomes relevan t for shock 
velocities >100 km s -1 dMicelotta et al.l l20H)al ). It is un¬ 
likely that the low velocity shocks that heat the molecular 
gas and enhace the H2 emission are also destroying the PAH 
grains, althoug h we clearly see a cha nge in the grain popula¬ 
tion. Similarlv. lContursi et alJ (I2013T) find that shocks do not 
dominate the [CII] and [OI] emission in the wind. Overall, 
the enhancement in H 2 /PAH is likely caused by H 2 exci¬ 
tation by shocks, although we cannot rule out a relatively 
small effect due to PAH destruction. 

I 11 dense molecular clouds, where the PAH destruction 
by shocks is less efficient, cosmic rays are the most efficient 
source of destruction of PAH molec ules, most ly because o f 
their bombardment by energetic ions dMicelotta et al.ll201llf . 


This is likely to play an important role in the M82 wind, be¬ 
cause the intense star formation in M82 is likely to enhance 
the cosmic r ay intensity by a facto r of « 5 compared to the 
Milky Way dRobitaille et al.ll2010l 'l. Theoretical studies sug¬ 
gest that CRs can destroy PAHs in dense molecular clouds, 
but in these environments the heating will al so be very ineffi¬ 
cient, causing a reduction of PAH emission. Mic elotta et aid 
d 201 l|j showed that only PAHs with a number of carbon 
atoms greater than « 250 have a survival lifetime greater 
than the circulation time scale of the outflow (~ 2 x 10 s yr). 
Thus, there may be a size-selection effect of PAHs in the 
outflow, which is roughly consistent with the size inferred 
from the band ratios (we mostly see PAH 6 .2/7.7 < 0.3, 
meaning an average size greater than 250 carbon atoms, see 
Figure [5J, and cosmic rays are likely to play a role. 

The calculations of Micelotta et al. 2010 show that any 
PAH shocked at > 100 km/s should be destroyed. To be 
re-formed in the post-shock gas, an efficient and fast re¬ 
formation process is needed. PAHs are a product of high 
temperature chemistry involving abundant carbon bearing 
precursors such as CH 4 and C 2 H 2 . In the cooled postshock 
gas these precursor species are likely not very abundant. 
In the cold neutral medium, the global PAH coagulation 
timescale is of the order of 100 Myr , muc h longer than the 
outfl ow dynamical time dSeok et al] l2014h . AKARI results 
from lYamagishi et al.l d 2012 l f suggest that very small grains 
still exist in the halo. Therefore it is reasonable to assume 
that PAH are produced in the disk, carried up in the outflow, 
and ’’protected” from fast shocks in the denser material. 
If molecular clouds are formed in the postshock gas in the 
outflow, then there may be PAH coagulation over timescales 
comparable to the outflow time (coagulation timescales are 
of the order of a few Myr in molecular clouds), but it should 
be very small fraction of the volume, contrary to what is 
observed. 

In conclusion, the enhancement of the H 2 /PAH is likely 
caused by H 2 shock excitation, although we cannot rule out 
a relatively small effect due to PAH destruction. 


5 CONCLUSION 

• We have mapped the Northern and Southern superwind 
in the nearby starburst galaxy M82 using the Spitzer-TRS. 
We find strong H 2 , PAH and Ne line emission that extends 
over nearly 5 arcminutes above and below the plane of the 
galaxy, and have measured the variation of the PAH, 
[NeIII]/[NeII] and PAH feature ratios throughout the wind. 

• We estimate a total warm (100 — 500 K) molecular gas 
mass to be between ~ 5 — 17 x 10 6 Mq and the total kinetic 
energy of the warm molecular outflow to be between « 6 — 
20 x 10 53 erg assuming an outflow velocity v 0 utfiow= 160 km 
s _1 . The warm molecular gas provides about 0.5-5% of the 
kinetic energy - comparable to that observed in the atomic 
gas. 

• There are clear variations in the PAH band ratios 
throughout the wind. The 6 .2/7.7 and 11.3/7.7 micron PAH 
band flux ratios, when compared to models, suggest a higher 
fraction of larger and/or more ionized PAHs in the wind 
compared to the starburst disk, especially in the Northern 
part of the wind. 

• The warm molecular gas to PAH ratio (H 2 /PAH) is 
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enhanced in the outflow by factors of 10-100 as compared to 
the starburst disk. This enhancement in the H 2 /PAH ratio 
does not seem to follow the trends with the ionization of 
the atomic gas (as measured with the [Nelli] / [Nell] line flux 
ratio) seen in the starburst. Our results indicate that UV and 
X-ray emission are probably not the dominant excitation 
mechanisms of H 2 . Models suggest the observed values could 
be reproduced with slow C- and J-shocks (v < 40 km s -1 ) 
driven into dense molecular gas entrained in the outflow. 
Using a dynamical timescale of the hU-rich outflow of about 
10 7 years results in an upper limit of 5x 10 9 Mq of molecular 
gas heated and shocked over the lifetime of the wind. Given 
the rate at which the gas must be shocked, the reservoir for 
this shocked gas is the cold molecular gas as observed in CO. 

• The emission from all of these phases is likely due to 
clouds entrained in and accelerated by the flow. They are be¬ 
ing shock heated by the hot plasma with high ram pressure 
generated by the collective thermalization of stellar winds 
and supernova in the starburst region of M82. These clouds 
are multiphase and the variation in properties are likely due 
to a range of shock velocities which depend on the density 
of the pre-shocked gas which depend on the density of the 
pre-shocked gas. 
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Table 1. PAH band and ionic line fluxes (in units of 10 16 Wm 2 ) and ratios. The regions for the extracted 
spectra are shown in Figure |3] and described in the text. 


Region 

6.2/im 

7.7/im 

8.6^m 

11.3/im 

17 fim. 

PAH 6.2 - 12.6/xm 

6.2/7.7 

11.3/7.7 

[Nell] 

[Neill] 

i 

7.12 

23.31 

4.19 

6.44 

3.52 

44.58 

0.305 

0.151 

3.62 

1.36 

2 

1.86 

3.31 

1.28 

2.56 

0.920 

9.93 

0.562 

0.278 

2.60 

0.966 

3 

1.01 

2.58 

0.526 

1.12 

0.751 

5.987 

0.391 

0.291 

1.27 

0.550 

4 

0.811 

4.76 

0.895 

2.65 

1.76 

10.876 

0.170 

0.370 

0.536 

0.893 

5 

12.53 

40.41 

8.93 

13.86 

8.32 

84.05 

0.310 

0.206 

2.17 

5.38 

6 

7.30 

22.03 

6.04 

10.65 

7.94 

53.96 

0.331 

0.360 

2.79 

4.18 

7 

9.39 

29.26 

6.61 

11.43 

7.47 

64.16 

0.321 

0.255 

6.06 

3.44 

8 

349.6 

1167.7 

254.3 

397.7 

282.6 

2451.9 

0.299 

0.242 

339.3 

89.4 

9 

741.9 

2783.7 

513.5 

668.9 

254.3 

4962.3 

0.261 

0.091 

568.0 

83.1 

10 

10.3 

37.6 

82.7 

18.3 

10.7 

159.6 

0.274 

0.285 

5.43 

2.57 

11 

30.9 

105.7 

22.2 

37.0 

21.3 

217.1 

0.292 

0.202 

30.2 

9.91 

12 

8.22 

27.2 

7.55 

15.0 

9.40 

67.37 

0.302 

0.346 

9.37 

4.17 
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Table 2. H 2 line fluxes (in units of 10 17 Wm 2 ) and ratios. ”Total” denotes the sum of all emission in 
the maps. 


Region 

H 2 S(0) 

H 2 S(l) 

H 2 S(2) 

H 2 S(3) 

H 2 S(l) - S(3) 

1 

< 0.313 

4.13±0.12 

0.8404=0.142 

1.704=0.14 

6.67 

2 

- 

4.354=0.15 

0.7844=0.173 

3.944=0.10 

9.07 

3 

< 0.169 

2.814=0.07 

0.5174=0.267 

1.684=0.23 

5.01 

4 

< 1.45 

6.684=0.14 

1.744=0.16 

3.174=0.14 

11.6 

5 

< 0.991 

14.34=0.2 

3.934=0.24 

8.744=0.41 

27.0 

6 

< 2.81 

9.684=0.26 

6.244=0.55 

2.054=0.18 

18.0 

7 

< 1.25 

8.874=0.14 

2.374=0.14 

4.794=0.20 

16.0 

8 

< 48.7 

37.84=0.6 

32.24=0.4 

83.44=0.7 

153.4 

9 

- 

179.84=0.8 

84.64=0.3 

108.44=0.3 

372.8 

10 

< 0.528 

23.94=0.2 

5.734=0.11 

119.04=3.0 

148.6 

11 

< 5.51 

31.74=0.2 

10.44=0.1 

16.64=0.1 

58.7 

12 

< 0.613 

8.854=0.12 

2.504=0.12 

4.944=0.21 

16.3 

Total 

< 253.9 

1171.54=3.2 

455.14=2.7 

1010.04=5.6 

2636.6 


Table 3. H 2 Temperature Model Parameters. ”Total” denotes the sum of all selected regions. 


Regions 

T(K) 

Ortho/Para 
xlO 19 cm -2 

N(H 2 ) 

x10 4 Mq 

M(H 2 ) 

xlO 4 Lq 

L(H 2 ) 

Region 1 

210. 4= 18. 

2.883 

2.12 

4.284=1.47 

1.43 


1094. 4= 208. 

3.000 

0.05 

0.01074=0.0036 

1.83 

Region3 

2404=54 

2.945 

0.838 

12.94=10.1 

7.91 


8194=126 

3.000 

0.007 

0.1134=0.087 

7.90 

Region4 

236. 4= 14. 

2.940 

2.07 

4.194= 0.82 

2.26 


831. 4= 124. 

3.000 

0.0128 

0.034=0.01 

1.84 

Region5 

206. 4= 17. 

2.8871 

6.77 

13.74=4.1 

4.32 


549. 4= 54. 

3.000 

0.132 

0.274=0.11 

4.37 

Region6 

166. 4= 21. 

2.653 

13.5 

27.34=18.2 

3.25 


593. 4= 56. 

3.000 

0.086 

0.1524=0.052 

3.30 

Region7 

224. 4= 12. 

2.917 

3.30 

6.664=1.17 

2.96 


657. 4= 75. 

3.000 

0.0382 

0.07714=0.0297 

2.68 

Region8 

324.4:12. 

2.994 

5.32 

10.74=0.6 

64.0 


595.4=13. 

3.000 

0.80 

1.64=0.2 

3.6 

Region 10 

218.4=6. 

2.913 

9.5 

19.24=1.7 

33.0 


637.4=62. 

3.000 

0.091 

0.184=0.06 

5.2 

Region 11 

233. 4= 3. 

2.934 

11.4 

23.04=0.81 

11.8 


623. 4: 14. 

3.000 

0.02 

0.394=0.03 

10.1 

Region 12 

241. 4: 10. 

2.946 

2.61 

5.284=0.666 

3.17 


816. 4= 157. 

3.000 

0.02 

0.042 4=0.002 

3.09 

Total 

2724=2 

2.976 

1.39 

4764=7.66 

463 


8544=18 

3.000 

0.023 

7.674=0.42 

640 
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Table 4. H 2 Shock Model Parameters. ”Total” denotes the sum of all selected regions. 


Regions 

Vshock 

km s —1 

Mflow 
Mo/ yr 

^cool 

xlO 3 yr 

M(H 2 ) 
xlO 4 Mq 

Type 

Region 1 

5.5 

42±3.1 

1.6 

6.84=0.5 

C 


35.0 

16±1.4 

3.4 

5.44=0.47 

J 

Region3 

5.5 

210±35 

1.6 

344=5.6 

C 


33.0 

150±11 

3.2 

484=3.5 

J 

Region4 

6.0 

48±6.1 

1.6 

7.64=0.96 

C 


35.0 

16±2.4 

3.4 

5.54=0.83 

J 

Region5 

5.0 

200±12 

1.6 

334=2.0 

C 


20.0 

1.9±0.14 

2.1 

0.414=0.03 

C 

Region6 

4.0 

1200±200 

1.7 

2004=33 

C 


21.5 

1.3±0.31 

2.1 

0.274=0.06 

C 

Region 7 

5.0 

130±32 

1.6 

214=5.2 

C 


21.5 

0.974=0.14 

2.1 

0.24=0.029 

C 

Region8 

5.5 

400±91 

1.6 

644=15 

C 


21.5 

174=4.4 

2.1 

3.64=0.93 

C 

Region 10 

5.5 

210±40 

1.6 

334=6.3 

C 


21.5 

1.5±0.24 

2.1 

0.314=0.05 

C 

Region 11 

5.0 

500±63 

1.6 

804=10 

C 


18.0 

6.2±0.26 

2.1 

1.34=0.05 

C 

Region 12 

5.5 

844=6.3 

1.6 

134=0.98 

C 


22.0 

33±1.4 

2.2 

7.54=0.32 

J 

Total 

5.5 

100004=2100 

1.6 

17004=360 

C 


21.5 

180=1=14 

2.1 

394=3.1 

C 



